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Abstract 

It was recently argued in jTj that black hole complementarity strains the basic rules 
of quantum information theory, such as monogamy of entanglement. Motivated by this 
argument, we develop a practical framework for describing black hole evaporation via uni- 
tary time evolution, based on a holographic perspective in which all black hole degrees of 
freedom live on the stretched horizon. We model the horizon as a unitary quantum system 
with finite entropy, and do not postulate that the horizon geometry is smooth. We then show 
that, with mild assumptions, one can reconstruct local effective field theory observables that 
probe the black hole interior, and relative to which the state near the horizon looks like a 
local Minkowski vacuum. The reconstruction makes use of the formalism of quantum error 
correcting codes, and works for black hole states whose entanglement entropy does not yet 
saturate the Bekenstein-Hawking bound. Our general framework clarifies the black hole 
final state proposal, and allows a quantitative study of the transition into the "firewall" 
regime of maximally mixed black hole states. 
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1 Introduction 



The question of how black holes absorb, store, share and release quantum information 
remains among the most vexing mysteries in theoretical physics [21 E]. Logical consistency 
of quantum black hole physics with the Bekenstein-Hawking entropy [H [2] imposes tight 
constraints on the microscopic theory. The discovery of the holographic principle [5] [6] and 
its realization in string theory [7] suggests that black holes satisfy all rules of quantum 
mechanics. The principle of black hole complementarity [8] is designed to reconcile the 
seemingly conflicting observations of an outside observer and an infalling observer. It 
postulates that 

i) black hole formation and evaporation are described via unitary quantum evolution 

ii) the region outside the stretched horizon is well described by QFT in curved space 

iii) to an outside observer, the stretched horizon behaves like a quantum mechanical 
membrane with microscopic entropy bounded by the Bekenstein-Hawking formula 

iv) an infalling observer can cross the event horizon without encountering any trouble. 

In a recent paper [T] , Almheiri et al formulated an interesting counter-argument against 
the simultaneous validity of these four assumptions. The reasoning of [1] , referred to as the 
"firewall argument", goes as follows. At late times the quantum state of the black hole is 
maximally entangled with the early radiation. Because entanglement can not be shared, 
and the information eventually has to come out, the late radiation can not be entangled 
with the black hole. Since the latter property is believed to be a necessary requirement for 
the quantum state to look like the local Minkowski vacuum near the horizon, one seems 
forced to conclude that the infalling observer can not pass through the horizon without 
experiencing dramatic events. This provocative conclusion has led to a lively debate |9J[10J. 

The firewall argument is the old information paradox in reverse. The original paradox 
was that locality and the smoothness of the horizon seems to imply that quantum mechanics 
falls short in describing black hole formation and evaporation. So if instead one adopts the 
rules of quantum mechanics and quantum information theory, smoothness of the black hole 
horizon becomes the real mystery pj] . The value added by the reasoning of PQ, and of other 
recent [121 EH EES EH] an d less recent [151 EE] works, is that it aims to bring the paradox more 
sharply into focus, by adopting a more precise quantum information theoretic language. 

The ramifications of the firewall argument are somewhat disturbing, since it could en- 
danger the validity of Hawkings derivation of the evaporation process, which stands at the 
very basis of this scientific discussion. Therefore, before accepting its conclusion, one should 
carefully examine the logic and assumptions that go into it. This requires a precise enough 
framework, that preserves all the essential physics and elements of the paradox. 
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Motivated by this challenge, we develop a systematic description of the evaporation 
process that keeps the rules of quantum mechanics manifest. The idea is to treat the black 
hole as an ordinary quantum system that interacts with a radiation field through absorption 
and emission, much like an atom would. We will not adopt postulate iv), that the horizon is 
smooth, but instead we will investigate if we can derive postulate iv) from the other three 
postulates. The argument of |T] indicates that this derivation may be possible for young 
black holes, that do not yet saturate the Bekenstein-Hawking bound, but should break 
down for old black holes that do saturate the bound. 

Besides the postulates i), ii) and iii), our approach relies on one additional physical 
assumption: we assume that the quantum mechanical amplitudes, that describe the transi- 
tion between different black hole states as a result of the emission of Hawking quanta, are 
ergodic matrices [17] . Their detailed form depends on unknown Planck scale physics, but 
their coarse grained properties are determined via standard thermodynamic reasoning. 

As another helpful tool, we will make use of the parallel, already exploited in [IjZJ, be- 
tween the physics of black hole radiation and the loss of coherence of a quantum computer, 
due to its interaction with a noisy environment. This analogy allows one to apply inventive 
theoretical techniques of quantum information theory [22l [23] , such as quantum error cor- 
recting codes (QECC's), to organize the internal state of the black hole [12]. QECC's are 
designed to counteract the quantum information loss. We will show that in the black hole 
context, the same technique can be used to safeguard the smoothness of the horizon. It 
enables the construction of microscopic observables, that play the role of Hawking partners 
of the outside radiation, relative to which the black hole state looks like a local horizon 
vacuum state. Smoothness of the horizon and the Hawking evaporation process are thus 
consistent with each other, and with the laws of quantum mechanics - at least for black 
hole states that do not saturate the BH entropy. 

QEC codes are not full proof [23]. As long as the black hole state has a less-than- 
maximal entropy, the code works with high fidelity: the local horizon geometry is smooth, 
and effective QFT is valid on both sides. When the black hole starts saturating the BH 
entropy bound, however, the QEC procedure becomes less reliable and sometimes fails - 
this are the instances when quantum information is leaking out [T2] . Since the code is also 
responsible for the reconstruction of semi-classical interior observables, the effective QFT 
description of the black hole interior breaks down in this limit. 

In section 2, we state our man assumptions and develop an interaction picture for black 
hole evaporation. In section 3, we present the reconstruction of the black hole interior via 
QEC technology, and analyze the firewall limit. For the well-informed reader that prefers 
the short version of our story, we have summarized our main arguments in section 4. 
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2 Setting the Stage 



In this section and the next, we will describe the basic setup and main assumptions. Our 
guiding principle is that a black hole is an ordinary quantum mechanical system, - much 
like an atom, albeit a very large one and with an enormous level density of states - that 
evolves and interacts with its environment via standard hamiltonian evolution. We start 
by setting up some notation that will help us model the evaporation process in a way that 
keeps the rules of quantum mechanics manifestly intact, in particular it obeys unitarity. 
Along the way, we present a version of the firewal argument, that can be directly addressed 
and tested within our framework. 

2.1 Black Holes States: Old and Young 

Consider a macrosopic black hole. Suppose that we have measured its mass M with some 
high accuracy. For now, we will ignore the uncertainty in M, and any other macroscopic 
properties such as angular momentum and charge. We associate a Hilbert space H M to the 
black hole with mass M, whose states 

\i) e H M (1) 

are indistinguishable from a macroscopic perspective. Each state | i ) describes the interior 
of a black hole, with microscopic properties that are hidden from the outside observer by the 
stretched horizon. The number of black holes micro-states is determined by the Bekenstein 
Hawking entropj^] 

N = dim H A / = e SBH , (2) 

with S BH = AnM 2 . We will assume that the time evolution of every micro state describes 
a black hole that slowly evaporates via the Hawking process. So as time progresses, every 
state | i ) G T-Lm will evolve into a new state given by the product of a micro-state of a black 
hole with smaller mass M' = M — E and a radiation state with energy E. 

To describe this evaporation process we need to enlarge our Hilbert space, so that it 
incorporates all possible intermediate situations. So we introduce the total Hilbert space 

■""We assume here that we have performed a suitable coarse graining, since strictly speaking the number 
of quantum states of a black hole should be represented in terms a density of states. This technical point 
will not be important for our subsequent discussion. 
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"H, given by the tensor product 



U = U BH ® H R , (3) 
where T-L BH is the space spanned by all black hole states with mass less than M 

Hbh = © H-m-Ei (4) 



E<M 



and 1-Lr is the Hilbert space of the radiation field outside of the stretched horizon. In the 
next subsection, we will set up a general description of the evaporation process as a unitary 
Hamiltonian evolution on the total Hilbert space (|3). Note that the black hole states | i ) 
are not part of the asymptotic past and future Hilbert space, but are to be regarded as 
very long lived resonances. It is standard in quantum mechanics to include resonances as 
an additional Hilbert space sector, even if they do not appear as true stable asymptotic 
states. Unitarity is preserved under time evolution, because the resonant black hole states 
eventually all evolve into the unique zero mass ground state |0) times an asymptotic 
radiation state. We first summarize what one would expect the time evolution to look like 
from a somewhat coarse grained perspective. 

First consider a young black hole of mass M. Let us imagine that it was formed in a 
pure state I i ). The total state of the black hole with its environment then takes the form 

|*young) = \i) |$>. (5) 

Here |$) denotes the state of the environment. Since the black hole is assumed to be 
formed in a pure state, there is no entanglement between the black hole and its environment. 
Tracing over the environment yields the pure state density matrix 

PIT" = \^)(^\■ (6) 

Next consider an old black hole of mass M that has formed in the far past and has evap- 
orated for a very long time. At this late stage, the black hole and radiation are described 
by a very-close-to-maximally entangled quantum state of the form 

1**0 = -7X7 (7) 



'N 

Here the sum runs over all basis states \i ) of the Hilbert space Hm- The precise form of 
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the states G "Hr depends on the early history of the black hole. Since the dimension 
of Hr is much larger than that of %m, the states |$«) can be safely assumed to form an 
orthonormal set. Hence when we trace the density matrix p = |^V^| over the radiation 
Hilbert space Hr, the resulting black hole density matrix p BH becomes, to a very good 
approximation, proportional to the identity matrix on the Hilbert space %m 

Pbh = ^EI*>M= ^m- (8) 

i 

In the following, we will consider the quantum states of both old as well as young black holes 
and take them as our initial condition at some time t — 0. Since hamiltonian evolution is 
a linear process one can express the density matrix of an entangled black hole at time t in 
terms of the evolved quantum states of initially unentangled black holes. We will make use 
of this fact to study the evolution of the density matrix p B H as the black hole continues to 
evaporate, and to see if it is possible to uncover an effective description of the semi-classical 
vacuum state near the horizon, as seen by an infalling observer. Note that even for young 
black holes this is not automatic, since we have not assumed the existence of an interior 
part. One of the main results of our paper is that this can indeed be done, with only a 
minimal set of assumptions. 

The maximally entangled old black holes and completely unentangled young black holes 
are the extreme limits of a broad class of intermediate cases going from slightly to almost 
maximally entangled situations. We will also consider these intermediate cases, since they 
are more realistic. Indeed one could argue that even young black holes can never be close 
to a pure state. As emphasized e.g. in [IB], in order for two different space time regions to 
be connected parts of the same space time, these regions have in terms of the microscopic 
physics to be maximally entangled. So in this sense, even young black holes may already 
come close to saturating the BH entropy bound. 



2.2 The Firewall Argument 

We now briefly summarize one version of the argument in pp. Consider the old black hole. 
Following pQ, we factorize the total Hilbert space into an early and late part 

U = U l&tc ®U™ ly . (9) 

Here %^ arly contains all radiation emitted during the early life of the black hole before t = 0. 
The late Hilbert space % late decomposes into the product of the Hilbert space of black hole 
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states with mass M' < M, times the Hilbert space of all radiation that will be emitted at 
late times t > 0. For notational convenience we will denote the two tensor factors of % late 
by 

y]ate _ ^ ^ __ % BH ^ 1{ B = %^ te . (10) 

Indeed, the firewall argument involves an Alice and Bob gedanken experiment. Alice passes 
through the horizon at some time t = r > and observes the interior black hole state in 
H A . Bob, on the other hand, initially stays outside and then falls through the horizon a 
small time later. He then catches up with Alice to confirm that she safely fell through the 
horizon. Suppose he succeeds and concludes that both are unharmed, for now. Why is this 
a problem? To explain the issue, we need a bit more notation. 

At t — t, the moment that Alice falls in, the black hole state ^ has evolved into a 
mixed state on the tensor product Hilbert space % late and also contains some radiation. 
Let us denote the orthonormal basis of late radiation states by 

| n) e H B . 

The states | n ) span the Hilbert space of a quantum field theory in the curved space-time 
background in the neighborhood of the black hole. In what follows, we will not need any 
detailed information about this QFT: all our arguments and calculations will go through 
regardless of whether it is strongly interacting or free. But for concreteness, consider the 
case of a free field theory. We can then expand the fields as usual in terms of creation and 
annihilation modes via (p(t) = ^= [b^e %ljjt + fr^e - ^*) . Our short hand notation | n) is 
then short-hand for the usual particle number eigenstates 

\^ = U^(>>ir\ o ), (») 

In particular, the states | n) 6 are energy eigenstates of the Schwarzschild Hamiltonian H. 

The central player in the firewall argument, and our main object of study, is the density 
matrix of the black hole and late radiation at time t = r > 0. It is defined on the tensor 
product Hilbert space H A (g) 7i B and takes the general form 

PAB = J2pnm\ n ) b ( m \- ( 12 ) 

n,m 
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The partial density matrix p nm acts only on Ti A , and specifies the correlation between the 
black hole state and external radiation. We do not have much a priori information about 
its detailed form, although in the following we will find out quite a bit about its coarse 
grained properties. Indeed, we have one important piece of knowledge: thanks to Hawking, 
or more directly (given that we already postulated the size of the BH Hilbert space) by 
application of standard rules of statistical physics, we know that if we take the partial trace 
of over the internal Hilbert space, p B = tr A (p AB ), we obtain a thermal density matrix 

PB = 2^w n \n) b (n\ ■ Wn =—^-, (13) 

n 

with Z = Y Jn e~^ En and (3 = 8ttM the inverse Hawking temperature. Here we have normal- 
ized the Boltzman weights w n such that 

^b(pb) = Xy„=l. (14) 



Note that equation (13), although it looks static, in fact represents a time dependent state, 



by virtue of the fact that the level density of allowed frequencies u increases inversely 



proportional with the time interval between and r. So (13) describes a radiation state 
with total average energy E = E n w n that increases linearly with time r. In the following 
we will assume that the time r is short compared to the black hole life time, so that E -C M. 



The statistical mechanical derivation of the density matrix (13) makes use of the fact 



that the number N n of microscopic black hole states after the emitting radiation equals 

dimU M - En = N n = Ne~^, 
while the total number of states in % A <g> % B with total mass M is given by 

dim (U A ® U B ) M = dimn M - En = NZ. (15) 



Equation (13) thus follows from the fact that each black hole state is occupied with the 
same uniform probability. 

Let us now return to the firewall argument pQ. The initial condition that the black hole 
starts at t = in the maximally entangled state ^ heavily restricts the possible form of 
the density matrix at the later time t — r of the black hole and the radiation that was 
emitted during that time. It was argued in [1] that, after letting the state evolve and emit 
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an extra segment of radiation, the maximally mixed form is essentially preserved. Suppose 
that this is true. One then arrives at the following form of p nm ' 



O _ — _ — ^ t / 16 x 

^nm M £ nm ^M-E n - \ ±KJ J 



Now comes the firewall argument: a density matrix of the factorized form (16) is clearly in- 
consistent with the assumption that Alice can fall through the horizon without disturbance. 
A direct way to see the problem is to compute the mutual information 

Jab = S a + S b ~ S ab ( 17 ) 



between the interior Hilbert space T-L A and the outside Hilbert space T-L B . Using eqns (|16|) 



and (13), one finds 



^ MPS = log(iVZ), S A = logAT- /3E, S B = \ogZ + 0E, (18) 

where E = Y2 n w nE n is the average amount of energy contained in the late radiation. Notice 
the appearance of the extra term log Z = —/3F, with F the free energy of the radiation]^] 
We will comment on its significance in a moment. 



Equation (18) tells us that the black hole interior and the late radiation share zero 



mutual information I^ lFS = S . + S D — S^ lPS = 0. A local vacuum state near the black 

AB A B AB 

hole horizon, on the other hand, encodes measurable entanglement between modes on each 



side. The density matrix (16) evidently does not. The authors of [lj argue that this 
discrepancy is a logical consequence of quantum information theory: since all black hole 
degrees of freedom are entangled with the early radiation, and entanglement can only be 
shared once, there is no room left for entanglement between the black hole interior and the 
late radiation. This contradiction led the authors of [lj to conclude that the horizon of an 
old black hole can not be a smooth region of space time. 

The above version of the firewall argument is not completely precise. In particular, the 



maximally mixed state (16) is not consistent with unitarity. To see this, note that the 



entropy = log NZ of the state (16J) is larger, by an extra term of logZ, than the 

BH entropy S BH = logiV of the initial state (|8]). The origin of the extra term is easily 
understood. The time evolution takes place in Hilbert space % late = Ha (g) which is 
much bigger than the interior Hilbert space Ha — T^bh by itself. The initial state ^ 

2 For a free photon gas, @F = —/3E/3, leading to S B = (4/3)(3E. The fact that the entropy in the 
radiation is a factor 4/3 bigger than the change AS^^f = —j3E in the BH entropy was first noted in [21] , 
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thus spreads out into an enlarged Hilbert space, and its coarse grained entropy gains an 
extra contribution proportional to the associated Helmholtz free energy F. Unitary time 
evolution, on the other hand, preserves microscopic entropy So (16) can not be exactly 



equal to the density matrix (|8j) evolved to a later time t = r. There must be off diagonal 
corrections that restore the unitarity. 

We thus have found a small gap in the reasoning of [1] . In [1] the black hole and emitted 
Hawking radiation combined were treated as if it forms a closed system, with a Hilbert space 
of some fixed size N. However, as the above simple analysis indicates, such a framework 
is too restrictive to describe the Hawking emission process. Eqns (13) and (12) only make 
sense as relations on a tensor product space % late = "H A <g> J{ B on which the b- modes can 
act independently. The local environment of the black hole, into which it slowly releases 
radiation and information, has a separate Hilbert space that is spanned by states that do 
not evolve from the black hole states themselves. A radiating black hole is not a closed 
system, but an open system, just like a radiating atom. Monogamy of entanglement holds 
in closed systems, but becomes less restrictive in open systems. 



2.3 Evaporation in the Interaction Picture 

Clearly, we need a more precise framework to model the evaporation of the black hole states 
(16]) and (l8j) in a way that manifestly preserves unitarity. We have already characterized 
a Hilbert space. We also need to characterize our Hamiltonian that describes the time 
evolution. We will choose to describe the time evolution in terms of the Schwarzschild time 
t, the time coordinate used by an outside observer. 

We first consider time evolution on the full Hilbert space, including the early radiation. 
In accordance with the decomposition ([3]), we write the Hamiltonian H as 

H = H BH + H R + H inU (19) 



where H BH acts only on interior black hole states, H R governs the evolution of the radiation 
outside the stretched horizon. H int represents the interaction Hamiltonian, that is respon- 
sible for the Hawking radiation process. The precise form of all terms in the Hamiltonian 
depend on unknown micro-physical details. We will therefore make only minimal assump- 
tions about the specific form of each term. All that we need to know is that H BH and H R 
each act within the respective Hilbert spaces, and that Hi nt acts on the tensor product in 
a way that allows the two sectors to reach a quasi-static thermal equilibrium. 

We can now use the familiar interaction picture, in which the time dependence of op- 
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Figure 1: A stationary element of the Hilbert space on which the interaction Hamiltonian acts 
is a tensor product of an eternal black hole state \i) A with a state representing the radiation 
outside the stretched horizon. The radiation state at t = factorizes into a product of a state 
I "Nearly con taining the early radiation sand the vacuum state | ) B of the late radiation. 

erators is governed by the free Hamiltonian H BH + H R , and states evolve through the 
interaction Hamiltonian H int . Following standard time-dependent perturbation theory, we 
may then expand the time-dependent states [^(t)) as linear combinations, with time de- 
pendent coefficients, of static eigenstates of the non-interacting Hamiltonian. The static 
states are of the form | j ) | n ) where \j ) G H BH is a microstate of a non-radiating eternal 
black hole of mass M — E n , that exist in the fiducial reference world in which the interaction 
Hamiltonian has been turned ofj^J and | n) G Hr represents a static radiation state with 
energy E n that lives outside the stretched horizon of the static eternal black hole. This 
sitation is depicted in the left diagram in fig. 1. 

Following p], we assume that the time evolution of states in Hr are well described 
by ordinary quantum field theory on curved space-time, and that all interactions with the 
black hole microstates take place at the stretched horizon. In the context of the interaction 
picture this means that the time evolution of the radiation field can be represented in the 
static black hole background. Thus, all complications of the evaporation process, including 
the backreaction, are contained in the time evolution of the quantum states generated by 
the interaction Hamiltonian Hi nt (t). 

To proceed, we need to take slightly better look at the radiation states Super- 
ficially, specifies initial conditions only for radiation on a spatial slice outside the 
stretched horizon. However, not all modes in the outside region need to pass through this 

3 In other words, \j ) is assumed to be an eigenstate of Hgjj with eigen value M — E n . 
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slice. As indicated in fig. 1, the stretched horizon is a semi-permeable membrane from 
which right-moving modes can emanate. Hence, the Cauchy data need to include initial 
conditions at or inside the future stretched horizon. The natural and only consistent choice 
is to impose vacuum boundary conditions, as measured in the Schwarzschild time, for the 
outgoing modes at t — 0. In this way we guarantee that the mass M — E n of the eternal 
black hole is all contained in the interior part | j ) of the wave function. 

From a semi-classical space-time perspective, it looks perhaps a bit suspect to assume 
that the radiation modes start out in their ground state as measured in the static frame, 
since this represents a singular state to an infalling observer. It is important to keep in 
mind, however, that smoothness of the horizon is not an input assumption of our set-up. 
Our basic principle is to preserve the rules of quantum mechanics, and to describe the 
evaporation process in the same way as one would describe the emission of a photon by 
an atom. With this starting point, Hawking radiation is not a consequence of a geometric 
regularity condition, but the result of turning on the interaction Hamiltonian Hi nt . From 
this perspective, the fact that the late radiation starts in the vacuum state means that one 
is dealing with spontaneous emission, rather than stimulated emission. 

With this motivation, we postulate that the quantum states of the radiation at t = 
factorize into the tensor product of state |$«) early , describing the early radiation, times a 
vacuum state | V that specifies the initial condition for the late radiation. 

NH°>J^W ( 2 °) 

Here | ) 6 denotes the Boulware vacuum of the late radiation. The state |$i) eariy is sta- 
tionary, since the early radiation is assumed to have already left the interaction region at 
t = 0. In other words, the early state is entirely decoupled from the late evolution. In 
the following we mostly forget about |^i) early - All we need to know is that after tracing 
p = j^V 1 !'! over "H^ arly we end up with a maximally mixed state of the form (8), multiplied 
with the vacuum density matrix |o) (0| for the late radiation. 

We are interested in the time evolution of the quantum state on "% late = % A <g> % B) 
describing the interior black hole together with the late radiation, starting from t = 0. It is 
sufficient to consider this time evolution for pure initial states | i )| ). This state evolves 
through the interaction hamiltonian H int , and thus after time r it takes the form 

U(r)\i)\0) b 
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where U(t) denotes the evolution operator 



U(t) = Texp (- % -JdtH int {t)}. (21) 

Since the interaction Hamiltonian couples the black hole interior to the radiation space Ht, 
the time evolved state decomposes as 

^)hO|o> 6 = E^|i>K- (22) 



n, J 



The coefficients C7 n l • represent the time-dependent probability amplitude for finding the 
black hole in the lower energy state combined with Hawking radiation in the state 



n)h, at the time t — r. The states | j ) in (22) are black holes with energy M — E n , which 



span a Hilbert space of dimension dim?{_M— e„ = e " n N. So for given n, the transition 
amplitudes C£j combine into a non-invertible matrix. 

In the following, it will be convenient to write the C n coefficients as explicit matrix 
elements C^j = (j\ C n \i) of operators acting on the black hole Hilbert space j^j 

C n = b (n\U(r)\0) b (23) 



Equation (22) can then be succinctly written as (here we simplify U{r) to U) 

U\i)\0) = J2 C n \i)\n). (24) 

n 

The unitarity constraint U ' U = 1 implies that 

J2 C ij C ti = E <%Cn = Uf. (25) 

n,j n 

Here and in what follows, C n and C\ will always denote maps from Hm^o Hm-e„ and back. 
For the rest, quite little is known about the j coefficients: their precise form depends on 
secret details of the micro-physics. Our main assumption will be that they look like large 
ergodic matrices |l7j . with statistical properties consistent with black hole thermodynamics. 
We will see that this will enable us to determine their coarse grained structure with sufficient 
accuracy to decode, with the help of some familiar quantum information technology |23j, 



'Note that, since U{t) acts on the tensor product T-Lj^ <8> 'H.Bi the rhs is still operator valued on T-La- 
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the local semi-classical physics on both sides of the black hole horizon. Equation (24) will 
form the starting point for much of the rest of our discussion. 

Our main story is quite independent of the detailed form of the interaction Hamiltonian 
H int (t). But for concreteness and to gain some intuition, it may be useful to have a specific 
model in mind. Suppose that the radiation is described by a free field 0. One of our 
postulates states that the black hole and radiation interact only at the stretched horizon. 
Locality then implies that Hi nt (t) is expressed in the radiation field <p(t), and its first 
derivative, at the stretched horizon. This leads to a simplified model of the form 

H int {t) = Jdt (<f>-Jt) V in (t) + 4> out (t) V out (t)) . (26) 

The interaction contains both incoming and outgoing fields, and hence describes black hole 
formation and evaporation. Here we focus on the evaporation process, which just involves 
the outgoing field (p out = <t>. The corresponding interaction operator V{t) can be expanded 
in modes V u and V^, that map black hole states with mass M to states with mass M ±tu. 
The matrix elements Vij(u) = (j V w |i) determine the emission and absorption rates of 
particles of frequency lo. Semi-classically this process is described by the creation of a 
Hawking pair. Hence it is natural to conjecture that the operators V(u) and V\uS) have 
something to do with the creating or annihilating particles behind the horizon. 

In principle one can express the C n operators in terms of the V operators as follows: one 
inserts the Fock states of the radiation and commutes the annihilation operators through 
the evolution operator. One obtains an expression involving products of a number of V 
operators, one for every emitted particle. It is suggestive that the number of V insertions 
matches the number of infalling partner modes of the Hawking pairs. This suggests that 
it is natural to view the transition amplitudes C n as creating a state with a collection of 
particles behind the horizon. Our aim is to make this intuition more precise. 

There exists of course a close parallel between the above discussion and the GKPW 
dictionary that underlies the AdS/CFT correspondence [7]. Indeed, we may compare the 
process of falling into a black hole with entering the near-horizon geometry of a stack of 
D-branes. In this comparison, the AdS space replaces the black hole interior geometry, the 
CFT represents the interior Hilbert space, the value of the scalar field at the stretched 
horizon plays the role of the asymptotic boundary value of the bulk field in AdS, and the 
V operators represent the CFT operator dual to 0. 
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2.4 Black Hole Entanglement Revisited 



As a first application of our set up, let us revisit the time evolution of the young and old 
black hole states (|6l) and 081). 



First consider the young black hole. At time r it has evolved into a state of the form 

Pr n§ = U\i)\0) b (0\(i\U\ (27) 



which we may expand in terms of the radiation basis as in equation (12), with 



p y :r = Cn\i)(i\ci (28) 

We now make the assumption that the amplitudes C n are large ergodic matrices. In 
practice this means that we can view their matrix elements as being randomly chosen from 
an appropriate ensemble. The properties of this ensemble are completely fixed by quantum 
statistics, and should be consistent with the usual laws of thermodynamics. 

In particular, if we start with a generic pure black hole state, the density matrix on the 



radiation states is expected to take the thermal form (13). Thus statistically one should 
have the following identity 



ClC n \i) = w n 5 mn) (29) 



with w n the normalized Boltzmann weight (13). This should be viewed as a coarse grained 
equality that holds for generic states \i ). 

Let us check this relation by writing the left hand side explicitly in components 

T,C i nij C^ j = w n 5 mni (30) 

3 

and use our knowledge about the size of the black hole Hilbert spaces. First take n = m: 
we are then summing Ne~^ En positive terms - so the answer is clearly non-zero - and since 
all amplitudes are statistically of equal size, the result must be proportional to e~^ En . The 



unitarity condition (25) tells us that the sum over n of the entire expression is equal to 1, 
so this fixes the normalization. Next suppose that n ^ m. Then we are summing of the 
order of Ne~@ En terms with random phases. This sum is strongly suppressed relative to 



the diagonal term. The relative suppression of the off-diagonal terms is of order 1/yN. 



We thus find that statistically the above identity (30) indeed holds 
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Before proceeding to the old black hole, let us determine the mutual information between 
the radiation and the black hole state. Since the initial state is pure, the entropy the black 
hole interior A, the radiation region B, and their union AB, is easily calculated, and gives 

5^ = 0, S B = \ogZ + f3E, 

young BH : (31) 

We see that the radiation region is maximally entangled with the internal black hole and 
carries the thermal entropy, as expected for a near horizon state. Indeed, on information 
theoretic grounds, there seems to be no obstruction to identify the young black hole state 



(27) with a smooth local vacuum state. We will make this identification explicit in the next 
section. 

Now let us look at the old black hole (|8l). After time r it has evolved into the state 



C = U P(°) U ^ with p(0) = ^lM®|0) fe <0| (32) 



The density matrix is obtained from that of the young black hole (27)by averaging over all 



initial black hole states. It can be expanded in the radiation basis as in ( 12 ) with 



Pnm — ~jy C"Cm- (33) 

The fact that the density matrix pb of the radiation region B is again thermal follows from 



tr A (C^C m ) = w n 5 mn which is automatically obtained from (29) by averaging over all 
black hole microstates. 

Let us again compute the entropies associated with the three regions A, B and AU B. 



First we note that the density matrix ( 32 ) manifestly satisfies 



(PB d H ) 2 = ^Pb% (34) 

from which we immediately see that it still carries entropy = log N, as it should. This 
result should be contrasted with the entropy S^ 1PS = log NZ of the diagonal mixed state 

Ad 



(16). The relation between the two density matrices can be understood if we introduce a 



phase averaging procedure, that averages over all possible phases of the transition ampli- 

n old 
AB ' 



tudes C- n . Performing this phase average eliminates all off diagonal terms of p°^, resulting 
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in the maximally mixed density matrix p 



AMPS 



1 



Pmn — C n C m — &nm^- M -E n - (35) 

The extra term log Z in the entropy is thus due to coarse graining as a result of taking the 
phase average. One should keep in mind that there are many off-diagonal corrections to 



(35). Individually, these are all small but there are many of them. So collectively, they can 



have a substantive effect. 



The diagonal form (35) is still appropriate for computing the entropy of the black hole 



interior region A. We thus find the following result for the entropies and mutual information 

S A = logN-PE, S AB = \ogN 
old BH : (36) 
S B = \ogZ + f3E, I AB =logZ. 

We see that the radiation region B and the interior region A now do share mutual infor- 
mation! This means that there is a non-zero amount of entanglement between the two. 

How is it possible that the maximally mixed state ^8J) evolves into a state with non- 
zero mutual information between region A and Bl Clearly, this would not be possible 
if we would think of Hawking particles as quantum bits that were previously contained 
inside the black hole. The key point is that the state ^ is embedded inside a larger 
Hilbert space that, besides the black hole states, also contains the radiation modes. Since 
monogamy of entanglement does not hold in open systems, the time evolution can generate 



new entanglement. Note however that, comparing with eqn (31), the amount of shared 
information is a bit less than for the young black hole. For a free photon gas, log Z = 0E/3, 
so the amount of entanglement for the old black hole is roughly a factor of 8 less than for 
the young black hole. The question we will need to investigate is if this amount of shared 
entropy is sufficient to describe an approximate vacuum state near the horizon. 

Above we have used quantum statistical arguments to determine sums of products of 
transition amplitudes C % n j and their complex conjugates. We will use this same reasoning 
repeatedly in the following sections. One may try to formalize this using a random matrix 
approach. One practical choice would be to pick a gaussian matrix ensemble, where one 
can use Wick's theorem. A more accurate characterization of the matrix ensemble is as 



the space of all matrices C l n ^ that satisfy the unitarity constraint (25). The language of 
Wick contractions still serves a useful terminology for labeling and estimatingthe different 
resonant contributions. We will not try develop this formalism further here, since for the 
cases of interest the relevant result can be determined by simple reasoning as above. 
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3 Passing through the Horizon 



In this section we will take a more detailed look at this question. Our goal is to first make 
explicit how a relatively young quantum black hole, given either by a pure initial state ^ 
or a partially entangled quantum state can in fact describe a smooth horizon state. From 
a quantum information perspective, this result is not a real surprise.. Nonetheless, it is 
still an important first step. It has up to now been far from obvious how any quantum 
mechanical model of a black hole can be consistent with having a smooth horizon. This 
logical tension is at the heart of the black hole information paradox. 

Specifically, we will show how the transition amplitudes C n can be used to reconstruct 
a mirror of the QFT Hilbert space, spanned by the Hawking pairs of the outside radiation 
modes. Our construction makes use of the close analogy between our description of the 
evaporation process and the loss of coherence of a quantum computer p21 E3] . In the case of 
a quantum computer, one can recover the lost quantum information via the use of quantum 
error correction (QEC). This new tool will then enable us to investigate the issues raised 
by the information paradox, and the firewall argument in more quantitative terms. 



3.1 Quantum Error Correction 

Many of the quantities and relations introduced above have direct cousins in the quantum 
information context. Let us briefly mention a few entries of the dictionary (221 [23] • 

One can think of the black hole as a quantum storage device exposed to environmental 



noise. Eqn (24) then describes the time evolution, the leakage of quantum information from 
the device into the environment, written in the so-called Kraus representation. The transi- 
tion matrices C n are known as Kraus operators. The emission of a Hawking quantum corre- 
sponds to the occurrence of an error, and the Kraus operators are typically denoted by E n . 
We will continue to use the notation C n . The evolution operator U = J2 n Cn ® l n ) b (0| 
is known as the error super operator, and the image of the C n operators is called an error 
subspace. 

Not all states in the Hilbert space can be simultaneously protected from error. A QEC 
code includes the specification of a code subspace, the space of encoded states that one 
wishes to protect. We will denote states in the code subspace by 

| i ) e n code . (37) 

The code subspace is chosen independently of the actual state of the system. In the black 
hole context, one may be tempted to try to identify the whole interior Hilbert space with 
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the code subspace, but this is not realistic. The presence of extra states and degrees of 
freedom outside the code space is essential, so that one can build in sufficient redundancy 
to safeguard the encoded information. 

We choose the encoded states to span a subspace within the Hilbert space 'Km of mass 
M, of dimension iV co d e <C N. We will comment later on the corrections that come into play 
when we try to take the limit iV CO de —> N. 

Not all types of quantum errors can be corrected. A necessary and sufficient condition 
is that for any pair of Kraus operators C n and C m , and for any pair of states | i ) and | k ) 
in the code subspace, one has 



C\C m \k) = w n 5 nm 5: 



k- 



(3* 



Here we fixed the normalization via eqn (29). If the property (38) does not hold, then 



errors can ruin the distinguishability of different encoded states and quantum information 



could be irrevocably damaged. For us, (38) is more than just a pre-condition. As we will 



see, it is the key algebraic identity that will enable us to decode the interior Hilbert state 
and exhibit the semi-classical horizon geometry. 



In the black hole context, the condition (38) is naturally satisfied: it generalizes eqn (29) 



to the case where |i) and /c) are different states, and follows from the same considerations 



outlined below eqn (29). It expresses the fact that C n and C m are statistically independent 



complex random matrices, that only interfere constructively when all terms in the compo- 



nent expansion of the product are of the form \C l n . ■] . Equation (38) should therefore not 



be read as an exact identity, but as an equation that for generic states holds with very 
high accuracy. Its corrections are suppressed by a relative factor of 1/yN. However, even 



while essentially true for any pair of states in H A , eqn (38) can not be read as a true global 
statement of the form C^C m = w n 5 nm , since this would appear to imply that C^C n is an 
invertible matrix, which it is not: its image has dimension N n = Ne~ l3En < N . So for given 



E n , we can use eqn (38) with confidence as long as we choose 



iV C ode N n = Ne 



We will encounter the necessity of this inequality at many other places. 

A quantum error correcting code is specified by the pair ("H C odc, {C n }), the code sub- 



space and the collection of errors it tries to correct. Given eqn (38), one can naturally 
associate a notion of entropy to a given QECC, given by Scode 



52 n w n Iogw n 
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Let us now describe the error correcting operation. It makes use of the recovery operators 



R n=^J2\ J )( J \ C n- ( 39 ) 

i 

These operators R n can be combined into a single super-operator R with the help of an 
'ancillary' Hilbert space T-L a , or simply 'the ancilla', spanned by basis states |n) in one-to- 
one correspondence with the basis states | n) b of the late radiation. Hence T-i a is isomorphic 
to H B . As an example, in case the radiation field is described by a free scalar field 0, the 
ancilla T-L a is the Fock space spanned by creation and annihilation operators a^ and a. 

The recovery super-operator R acts on the tensor product Ha <8> Ha via 

R \3)\ )a = J2 R n\j)\ n )a (40) 
n 

The ancillary Hilbert space is an important ingredient of quantum error correction, since as 
we will see shortly, it provides the depository into which it dumps the entanglement between 
the code space and the environment |23j. We emphasize, however, that in our context - 
unlike for the usual application of QEC codes - the ancillary Hilbert space 7i a is just an 
intermediate algebraic device, that helps us put the recovery operator in a convenient form. 
Since our main purpose is to keep track of information flow, it is therefore important keep 
in mind that the ancillary Hilbert space never really exists and needs to be projected out 
in the end, in a way that manifestly preserves unitarity. We will return to this point, once 
the physical meaning of the ancilla has become clear. 

The recovery operator is designed to reverse the error: applying R to a time evolved 
state U \ i )|0) a , one recovers the original state | i ) times some state on the Hilbert space 
T-i B ® r Ha, provided that the initial state | i ) lies inside the code subspace. The verification 



of this claim is straightforward, and makes use of the property (38) of the C„'s 
RU \^) a \H) b = J2 R m C n\i)\m) a \n) b 



^Y k )( l \ c l c A*)\ m ) a \ n \ ( 41 ) 



m,n,k 



5 Hcre we keep the short-hand notation U and R for the extensions U t a and R ® of the super- 
operators to the tensor product Ha ® H a ® %b- 
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We read off that 




if | % ) G Hcode 

RU\i)\0) a \0) b = { ' (42) 

if | l ) £ "Hcode 

where \ §u) is the following familiar looking state, defined on the tensor product of the 
Hilbert space of the outside radiation region with the ancillary Hilbert space 

\0u) = J2^~n\n) a \n) b . (43) 



Given that w n are the Boltzmann weights, we recognize this state as the Unruh vacuum 
state, with the ancilla states playing the role of the interior radiation region. We see explic- 
itly how the R operation works: it sweeps all entanglement between the time evolved black 
hole and the emitted radiation, encoded in E/|i)|o),, into entanglement of the radiation 
with the ancilla. Thus the recovery operation restores the purity of the black hole state, 
provided it started out within the code space. 



We see from (42 ) that the recovery operator R is not a unitary map, since it projects out 



all states that do not map back into the code space. This means that the operation R as 



defined in (40)-(39) in fact is not really a valid superoperator. Preferrably, one would want 



the recovery operation to be a unitary operation. To achieve this, consider the operator 

n = a <0| J R t J R|0) a . (44) 

We may also write II = ^2 n RnR n = Yli nW~ C n \i)(i \ Cn- If R had been a valid su- 
peroperator, then II would have been equal to the identity operator. Instead it defines a 
projection operator on the space of all states that can be reached by acting with the C n 
operators on the code subspace. In other words, it is the projection of all states that the 
code states can evolve into as a result of time evolution over time r. 



A straightforward calculation, similar to the one done in equation (41), shows that II 
indeed has the property of a projection operator: 

n 2 = n. (45) 

To complete the recovery super operator, one thus adds one extra term to the right-hand 
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side of (40) proportional to the projection operator 1 — II on the orthogonal complement. 




if | i ) E U coAe 

(l-n)l7|<>|0> 6 = { (46) 

if | % ) ^ H C ode 

We see that the QECC has indeed succeeded in its task of protecting all states within the 
code subspace from loss of coherence. 



3.2 Construction of the Interior Operators 

We have seen that the combined operation R U, time evolution followed by quantum error 
correction, associates to every black hole initial state \i ) G ^code an Unruh state of the 



form (43). Since for any state | i there is some recovery operator R for which | i ) lies 
within its code subspace, we have shown that | i ) (when evolved forward in time) always 
contains a factor that can be decoded into a local vacuum state near the horizon. 

We now define interior operators, that form the Hawking partner of the outgoing radia- 
tion modes. Specifically, we associate to every operator on H B a partner operator that acts 
on the interior black hole states. Our construction works for every state in the code sub- 
space. Suppose that the black hole together with the radiation at time t = are described 
by some mixed state of the form 

p(o) = j2 pik\ i )m ® i°>6<°i> ( 4? ) 

i ,k 

where the sum runs over some small subspace "Hsmaii of the total black hole Hilbert space. 
We are then free to choose a code subspace such that % S maii C T-L co de- At time t = r, the 



state (47) has evolved into the density matrix 

p(r) = Up(0) U ] . (4£ 



We wish to show that this density matrix describes a black hole with a smooth horizon. 

We will now define the interior operators. In the previous section, we introduced the 
ancillary Hilbert space 1~iai &s a carbon copy of the outside radiation Hilbert space H B . For 
us, this ancilla does not really exist, but is just introduced as a helpful algebraic tool. Even 
so, we are free to interpret % a as representing all states on the hidden Rindler wedge, the 
spatial region on the other side of the horizon from ~H B . Now consider some operator A 
that acts on 7-L a . We wish to identify A with the Hawking partner of an identical operator 
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acting on the outside radiation space 7i B . But we can not use the ancilla as our Hilbert 
space, so we need to transform A into an operator acting on the actual internal Hilbert 
space of the black hole. We now associate to A the following operator acting in the interior 
black hole Hilbert space 



a (0\R^AR\0) a . 



(49) 



where R is the recovery superoperator defined in (39)-(40). Note that this is a vacuum 



expectation value in % a . So A is a proper linear operator acting on 1-L A , the state space 



of the black hole, and the ancilla is indeed just a virtual device. Eqn (49) is modeled 



after the definition (44) of the projection operator onto the code subspace. Indeed, n 



assumes the role of the identity operator on the interior Hilbert space T-L a . 

As an example, in case the radiation field is described by a free scalar field <f), we can 
choose to consider A = <f)(y), where y is a point on the left Rindler wedge behind the 
horizon, that is, <p(y) is defined as a mode expansion in terms of ancillary oscillators and 



a. The general definition (49) then associates an effective local field operator 4>(y) } acting 
on the interior black hole Hilbert space, via 



4>(y) = a (0\R!<P(y)R\0) 



(50) 



Eqns (49)-(50) are our proposal for a microscopic realization of the interior QFT operators. 



To support this identification, we first need to verify that the map (49) between A 



and A is an isomorphism between operator product algebras. In other words, for any two 
operators Ai and A 2 , we need to show that 



A 1 A2 = a (0\R^A 1 A 2 R\0) 



(51) 



This equation tells us that applying the map (49) after taking the product gives the same 



result as taking the product after applying the map. If this holds, then the A operators 
indeed have the same operator product algebra as the effective QFT operators. 

Let us check that (51) indeed holds. f\ It is useful to expand out the definition (40) of 
the recovery superoperator. The definition ((491 can be expanded as 



(52) 



A = J2 R n A nmRr 



In this section we assume that iV co d <€. N. We will look at the case N co a c — > N in the next section. 
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where A nm = a (n\A\m} . We see that the combination R^Rm takes over the role of 
the operator |n) (m| in the ancilliary Hilbert space. As a free field example: with this 
notation, we may define the effective interior creation and annihilation operators a^ and a 
via 

a^ = Vji?^n|o^|m)i? ro and a = R\ l (n\a\rnjR m . 

n,m n,m 

The a and at operators act purely within the internal black hole Hilbert space T-i A , and 
thus commute with the creation and annihilation modes W and h of the radiation. 



With the definition (52), the requirement (51 ) that A satisfies the QFT operator product 



algebra amounts to the identity 

R) n R s R\R m = 5 st B) n R m - (53) 

As shown in the Appendix, this property indeed holds with great accuracy, as long as 
A^ C ode "C N. The precision with which it is valid is a measure for the successfulness of the 
quantum error correcting procedure R. For us, it measures the faithfulness with which we 
can recover the operator algebra of the interior low energy effective field theory. It is a 
strength of our formalism that we can in principle compute the corrections and limits of 
validity of the effective theory; we will consider the correction to the operator algebra when 
discussing the limit N code /N — > 1 in the next section. 



Next we compute the expectation value of the operators in the state p(r) given in (48). 



Here we assume that the initial state (47) at t = is contained inside the code subspace, 



and that A code <C N. Let A and B denote some general operators in terms of the a and 



6-modes, respectively, and let A be the interior operator defined via (49). The computation 
of the expectation value tr(p(r) AB) is straightforward, and summarized in the Appendix. 
We find that the expectation value coincides with that of the effective field theory in the 
local Unruh vacuum state 

tr(p(r)AB) = (0 V \AB\ V ). (54) 

This is the main result of our paper. It shows that the black hole state can look like a 
smooth horizon state, without straining the rules of quantum information theory. 



Another way to obtain this result is to use the expansion (52) in terms of matrix ele- 
ments. In the Appendix it is shown that 

tr A (p RlRm) = \/w n w m \ n) b (m\. (55) 
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Figure 2: Steps in the quantum computation of the expectation value tr(pAB) for a pure state 
density matrix p(0) = • The initial state |^)|0), first evolves via the evolution operator U . 
The interior operator then acts via a recovery operation R followed by the operator A acting on 
the ancilla. The outside operator B acts directly on the time evolved radiation state. At the 
instant that A and B act, the interior state is identical to the initial state and the a and b 
states are in the local vacuum state Oj/). Then follows the conjugate of the recovery operation 
a (o|i?^, which eliminates the ancilla. Finally, one projects back onto the final state t>(o|(z| . 



From this one directly computes (here B nm = b{n\B \ m) b ) 

tT(pAB) = Y, A nMp R] n R mB) = V^n^A nm B nm . (56) 



This is the expectation value in the local Unruh vacuum (43). 



What made this magic work? The definition (49 ) of the interior modes involves a QECC 
operation, which requires detailed knowledge of the emission amplitudes C n , whose precise 
form depends on unknown details of Planck scale physics. It looks like cheating to assume 
that Alice can access these observables, without actually having such detailed knowledge. 
The situation is of course similar to other holographic dualities, like AdS/CFT. From the 
dual perspective, the interior space-time emerges in a mysterious way. But for Alice, the 



smooth local space-time is a classical reality, relative to which the operators (49) behave 



like ordinary observables. So she has no trouble using them. Note, however, that Alice 
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does not organize the modes in terms of their Schwarschild energy u, but reassembles them 
into local Minkowski modes, that make no reference to the precise location of the horizon. 

To make this a bit more concrete, consider the two point function of local field operators, 



say one behind and one in front of the horizon. The general result (54) gives that 



tr(p(f) a (y)(j) b (x)) = (Ou\(f>a(y)M x )\ Q u), (57) 



where the subscript indicates that <\> a is made up from a-modes, and 4>b from 6-modes. 
What is the meaning of this expectation value? Strictly speaking, it does not yet describe 
what we want. In our discussion so far, we have only included the outgoing modes that 
are emitted by the black hole. The ancillary a modes play their guest role as the Hawking 
partner of the outgoing b modes: they initially live on the left region in fig 1, but can be 
analytically continued to live just behind the future horizon. To really get the complete 
set of modes behind the horizon, we would need to extend our framework a bit, by adding 
the infalling modes. Since the horizon state looks smooth, the infalling modes can be 
analytically continued to the region behind the horizon. So we have indeed explained how 
a QFT mode, and by extension any object, can fall into a black hole. 

Note that we can also try to go beyond the past horizon, by switching time direction. 
The outgoing modes behind the past horizon are obtained by analytic continuation of 
outgoing modes on the right region in fig 1. The time-flipped Hawking partners of the 
incoming modes live behind the past horizon, and are obtain by analytic continuation of 
modes on the left region in fig 1. To define them, we would need to introduce transition 
amplitudes for the ingoing modes involved in the black hole formation process. 



3.3 Ancilla versus Black Hole Final State 

Our use of the ancillary Hilbert space T-L a is somewhat reminiscent of the black hole final 
state proposal of [H]. However, there is a key difference: unlike the final state proposal, our 
procedure is completely in accord with the conventional rules of quantum mechanics. We 
have stated explicitly in what Hilbert space we are working, and time evolution is unitary. 



Our internal operators (49)] involve an expectation value inside an auxilary Hilbert space 



1-L a , but we could have avoided that altogether by directy using the definition (52). So we 
work at all times in the same Hilbert space, without tensoring in or projecting out any 
extra degrees of freedom. 

The black hole final state proposal of [19] works differently. There the ancilla T-L a and 
the radiation space Tib are new Hilbert space sectors that arise out of nowhere whenever a 
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back hole gets formed, and are postulated to start out in the state |0[/) = ^2 n V^n\ n } \ n }h- 

\i)^\i)\0u). (58) 

This is a crude prescription, however, that immediately obscures the rules of quantum 
mechanics. Where was the Hilbert space T-L a (g> Tib before the black hole was formed? How 
did it end up in a maximally entangled state \0u) with the radiation modes? How can one 
keep track of quantum information if such things can happen? 



Superficially, our equation (42) looks very similar to (58). For us, however, equation 



(42) does not represent time evolution, and the state | i )|0j/) does not exist in our physical 



Hilbert space H A £g> Tib- So let us interpret the proposal of [19] in the same spirit. In 



equation (58), we should not think of the ancilla states |n) as real, but virtual. Indeed we 
can undo the damage done by the rash introduction of the extra Hilbert space and restore 
unitarity by writing^] 

U\t)\0) b = a (0\tf\t)\0u). (59) 

The left-hand side is a manifestly unitary evolution of the initial state |^)|0) 6 , while the 
right-hand side looks like a non-unitary final state projection a <^0|i2^ applied to the state 
| i )|0[/). The apparent non-unitarity of the final state prescription cancels out against the 
initial unitarity violation of introduction of the additional Hilbert space sector. 

This interpretation immediately circumvents the main substantive objection against the 
final state proposal made in [20], where it was pointed out that any non-trivial dynamics 
between the two sectors would immediately destroy the coherence of the physically observ- 
able state. But as the above discussion makes clear, the ancillary Hilbert space should be 
viewed as dynamically completely decoupled from the physical Hilbert space: the evolu- 
tion operator U and Hamiltonian H acts only on 7-L A ® Tib and leaves H, a inert. That 
the dynamics exactly factorizes in this way would look very unreasonable if 7i a had been 
physically real, but is obvious in our set-up. The final state proposal is only salvageable - 
but then loses much of its explanatory potential - by realizing that the ancillary Hilbert 
space H a is an artificial construct. 

If the ancilla states are not real, how then are we are able to still employ quantum 
error correcting technology? QEC involves the reversal of decoherence. When the code 
block interacts with its environment, the two become entangled and the Von Neumann 



entropy of both systems increases, as seen from (31) The recovery operation R implements 



7 This equation follows from U\i)\0) b = TLU\i)\0) b = a (0\R f RU\i)\0) a \0) b = a (0\R f | i ) \0u)- 
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an entanglement swap: it purifies the code space by sweeping all its entanglement with 
the radiation modes into the ancilla. Hence reality of the ancilla is absolutely essential for 
quantum error correction, since where else can the entanglement entropy go? 



The resolution of this apparent conflict is as follows. As seen from the definition (49), 
-A = a (0|-R^4i2|0) a , the interior operators make only temporary use of the recovery proce- 
dure. Reading from right to left, they tensor in the auxiliary ground state |0) , and perform 
the recovery operation R. This step removes entropy from the code block and temporarily 
stores it inside T-i a . Then it applies A, and instantly undoes the recovery by acting with the 
'final state projection' a ^0|-R^. This final state projection - which is necessary for restoring 
unitarity - erases the ancilla and dumps the entropy back into the interior Hilbert space 
7-L a ^\ So the QEC procedure is able to perform an entanglement swap, but at a price: the 
interior black hole Hilbert space space clogs up with entropy as time passes by. 

This looks like a recipe for trouble. As the entropy of the black hole density matrix 
increases, one needs a bigger and bigger code space for performing the recovery procedure 
on the whole density matrix. Eventually N co d e starts approaching the size N of the full 
Hilbert space. But as we have seen, the quantum error correcting code becomes less and less 
reliable with increasing N co d e /N, and as a result, our identification of the interior effective 
QFT operators starts breaking down. This is a manifestation of the firewall problem. 



3.4 Fidelity and Errors 

Why do quantum error correcting codes have any useful application in the study of black 
hole information and entanglement? In the context of their application to quantum com- 
puters, QEC codes are designed to protect the quantum coherence of a particular part of 
the Hilbert space, the code subspace, by correcting the errors that occur by emission of 
particular quanta. It does so by undoing the transformation that is responsible for the 
emitted quanta of radiation, and moving the entanglement with these quanta to another 
part of the Hilbert space that is outside the code subspace. 

In the black hole context, the quantum recovery operation acts in an exactly identical 
way on the Hilbert states, but the QEC code plays a rather different role [12J. The code 
subspace is typically already in some highly entangled with the early radiation outside the 
black hole. The error correcting procedure still works by undoing the transformation due 
to the emission of radiation. If it succeeds, it entangles the radiation again with another 
part of the Hilbert space. As we have just shown, this entanglement helps preserve the 

8 This entropy transport goes hand in hand with transport of heat. In light of Jacobson's arguments 
[24j . relating entropy and heat flow with the GR equations, it would clearly be very interesting to study 
this in more detail. 
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entanglement across the smooth the horizon. 

In any QEC procedure there is always a small but finite chance that it fails. In the 
situation of quantum computer in an approximate pure state, the error will result into 
some partial decoherence: a small part of the code subspace will become entangled with 
the outside. As a result, the corresponding internal quantum state is no longer suitable for 
doing quantum computations. In the black hole application, failure of the QEC code means 
the opposite. Occurrence of an error means that the code did not preserve the entanglement 
between the code subspace with the early radiation. Therefore the QEC error will cause a 
swap in entanglement: it releases a qubit of information from the code subspace and lets it 
escape to the outside. This is how the information leaks out of the black hole [T2] . 

What is the fidelity and how big is the error of the QECC? To determine this, we 
need to look at how well the recovered state approximates the original state. A good way 
to compute this error is to determine how well the projection operator II, when evolved 
backwards in time, approximates the identity operator on the code subspace 

IH7|?)|0) b = U\i)\0) b + error. (60) 

This equality is without error for a perfect QECC with 100% fidelity. Plugging in the 
definitions for the left-hand side gives 

^ C n \k)(k\CtC m \i)\m) h . (61) 

k,n,m 

We evaluate this expression via the same procedure as before: we keep only those terms in 
which every matrix element ■ constructively interferes with its complex conjugate. In 
the limit iV co de <C N n , the dominant interference occurs between C^C m for the term with 



n = m. This leads to the first term on the right-hand side of (60). This is the good term. 

But there is also a bad term: when N co( x e /N n becomes non-negligible, the other inter- 
ference term becomes important: on the right can constructively interfere with C n on 
the left. In this way we find that the error is given by 

error = EU\i)\0) b , with E = £ ^ t M ^ En . (62) 

n n 

The operator E will be called the 'error operator'. Hence to have an accurate recovery 
procedure, we must choose the code subspace small enough so that this error is negligibly 
small. Thus we again encounter the condition that iV co dc <C N n . 
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3.5 Approaching the Firewall 



The firewall argument concerns the limit in which the black hole density matrix approx- 
imates the maximally mixed state ([8]). There are several motivations to study this limit. 
As emphasized in [TS], two space time regions can be connected if they are maximally en- 
tangled at a microscopic level. These arguments suggest that to be able to fall through the 
horizon without drama, one actually needs the maximal entanglement between the inside 
and the outside. Other related perspectives are the interpretation of black holes as informa- 
tion mirrors [T2Jand fast scramblers [T7], and Mathur's argument that, as far as quantum 
information goes, order one corrections are needed to the usual semi-classical picture. 

Motivated by these ideas, we will now study what happens as one approaches the max- 



imally mixed limit. We will still assume that the initial state p(0) in (47) fits inside a code 
subspace of size N code < N, which we parametrize via 



N, 



code 



e"^. (63) 



N 



The maximally mixed case corresponds to setting N co ,i c /N = 1, so that S co d e = logiV co dc 
saturates the Bekenstein-Hawking bound S BH = logiV. We will call this the firewall limit. 

Based on our set up, we can try to give a somewhat more precise characterization of 
the density matrix of an old black hole, as follows. As pointed out above, the quantum 
error correction procedure involved in reproducing the internal effective field theory is not 
perfect, but has a finite probability of producing an error. In the previous section, we saw 
that, when considering operators that act on an internal state with energy M — E n , this 
error is of the order 



N, 



code 



The occurrence of an error coincides with an event in which the black hole releases infor- 
mation about its interior state into the environment. So an important first question is: 
what is the expected maximal rate at which these errors can occur? At the late stages 
of evaporation, the number of information carrying quanta should exceed the number of 
quanta that increase the black hole entropy. Therefore, the error rate eventually exceeds 
the frequency of good recoveries. Our proposed definition of the density matrix of an old 
black hole is that its entropy has grown to the level such that the probability that any 
given Hawking particle carries out a qubit of quantum information into the environment is 
of order 1. This amounts to an error rate of about one qubit per black hole crossing time 
M. So in the following we will set our time step r equal to one crossing time. 



30 



The errors of the QECC give rise to corrections to the internal low energy effective field 



theory. Consider the operator combination W n R m featuring in the definition (52) of the 



internal operators. The erroneous Wick contraction produces an extra contribution 



RlR m = —== TC n \i)(t\Ci = 5 nm e l3 ^-^l M _ En . (65) 

i 

Here the bar indicates the phase average. The normalization of the right-hand side is 



determined by unitarity, c.f. eqn (35). Note that it indeed vanishes when e ^ — > 0, i.e. 



when the code space becomes negligibly small. This extra interference term deforms the 



algebra (53), that underlies the isomorphism between the internal operator algebras and 
that of the radiation field. Thus we need to reinvestigate our definition of the interior 
operators, and pick the definition that best preserves the effective QFT description. 

The most reasonable prescription is to introduce a normal ordered product :i?^i? m :, 
defined as the operator product with the phase averaged contribution RnR m subtracted. In 
other words, the normal ordered operators are traceless. This recipe preserves the identity 
tr A (p :R) m R n -)j = yjw n w m |n),(m|, which ensures that, if we keep our definition of 
the internal operators of the form (52) but now as normal ordered expression : A :, the 
expectation values of the effective QFT still look like those in the local vacuum state |0[/). 
So from the semi-classical perspective, it looks like the horizon is still smooth - in the sense 
that all expectation values are well behaved. 

The error has not gone away, however: the normal ordering procedure still leads to a 



deformation of the operator product algebra (53) 



:RlR s : :R\R m .= 5 st (:RlR m . + 5 nm e^ E ^H M . En ) . (66) 

To verify this relation, observe that the product of two normal ordered operators is not 
normal ordered; this leads to an extra contribution equal to RnR m . This extra term can 
be ignored as long as fi ^> E n . When the energy E starts approaching the value of fi, 
however, it starts to show up as a manifestation of the error of the QEC operation. Before 
we investigate the form of the deformation, let us determine how big or small it is. 

What is the probability that the QEC procedure produces an error? Consider the 



operator II introduced in (44). It can be expanded as II = J2 n : R^Rn'., and in the low 
energy effective theory, it represents the unit operator 1. It satisfies II 2 = II provided that 
the recovery precedure works perfectly. But once we include the possibility of errors, II is 
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no longer a perfect projection operator. Using (66), we find a correction term given by 



U 2 = U + E, E=J2^ {En ~^ )l M-E n . (67) 

n 

We assume that e _/3/i is adjusted so that E is still small. We may then interpret E as the 
operator that produces the error in the recovery operation. 

To determine the relative magnitude of the 'good' first term and the 'bad' error term, 
let us take their expectation value. For concreteness, we assume that the black hole density 
matrix at time t — takes the form p(0) = ^ - \ * )(* |- We then have 

tr A (p(r)U) = \ n ) b ( n \ = \z~ mb (68) 

n 

tr A (p(r)E) = Y,e-^\n) b (n\=e-^l B . (69) 

n 

The good term is the thermal density matrix p B of the radiation, and thus its trace is 
equal to 1. So it is normalized as a probability measure. When there is a non-zero error 
rate, however, this normalization needs to be adjusted, to account for the fact that the 
QEC operation is no longer has 100% fidelity. It seems reasonable to interpret the second 
equation as the relative probability that an error takes place in the quantum channel. It 
tells us that the errors have an infinite temperature: every state |n) in the radiation Hilbert 
space has , regardless of its total energy E n , equal probability p ~ e _/?M to be the carrier of 
some irretrievable error. This probability increases linearly with the iV co d e , as expected. 

As we increase the size of the code space, we see that there is a natural cross over 
point, where the error probability becomes of order 1. This crossover takes place when the 
dimension of the code space reaches the critical value 

Ncodc - P -P» ~ 1 (70) 
N ~ dimH B 1 } 

One may think of this as a detailed balance equation: it picks the saturation value for the 
size of the code subspace, such that the QEC procedure starts producing approximately 
one error per crossing time. Notice that the right-hand side is a UV sensitive quantity: it 
depends on the size of the short-distance cut-off of the effective field theory that describes 
the Hawking radiation)^] and on the location of the stretched horizon. 



3 We implicitly assume that the Hilbert space T~L B already comes with some natural IR cut off, say, on 
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When the error correction works, we get a semi-classical description of the evaporation 
process. But when it fails it produces modifications to the low energy effective field theory. 



What do they look like? Consider the identity (51 ) that verifies that the mapping (49)-(52) 
preserved operator product relations. We now find that the right-hand side receives an 
extra diagonal term, which we may think of as a result of normal ordering. We find that 
the operator product algebra is deformed according to 

:A V .:A 2 : = :A X A 2 : + ^(^^j^e^^M^, (71) 



How would a low energy observer interpret this deformation? It is tempting to look for 
a way to translate this algebra in pure effective field theory language, perhaps with some 



non-locality. However, the extra term is not of the form (49), and acts on a bigger Hilbert 
space than % co d e , the space we used for our re-construction of the effective QFT. Because 
of this, we have not found a sensible way to answer this question. 

Perhaps this is indeed the wrong question to ask. The presence of the extra term is 
dictated by unitarity of the underlying microscopic theory, and the microscopic mechanisms 
that are responsible for releasing the information from the black hole are likely to be not 
describable by semi-classical physics. In our QEC terminology, they produce error terms 
that in magnitude are of order one. Moreover, once the black hole entropy starts to decrease, 
the number of information carrying radiation modes must exceed the number of modes that 
increase the entanglement, the error rate must exceed the frequency of successful recovery 
operations. This indeed indicates that in this regime the semi-classical description can no 
longer be trusted. We are then led to conclude that when corrections to the semi-classical 
theory becomes important, the effective QFT disintegrates. This is a manifestation of the 
firewall problem, and perhaps the manifestation of the formation of an actual firewall or 
fuzzball. |25j . 

This last conclusion is still a bit premature, however. It is quite conceivable that the 
breakdown of the effective QFT is simply a consequence of our attempt to capture a close- 
to-maximally entangled state of a large system in terms of a single semi-classical reality. 
Indeed, what we have shown is that a black hole with small enough non Neumann entropy 
relative to its BH entropy bound, i.e. with e~^ M <C 1, supports a semi-classical geometry 
with a smooth horizon. Since a maximally mixed black hole is an incoherent sum of black 
hole states with less than maximal entropy, we have shown that it represents an incoherent 
sum of different (in the sense that they are incompatible with each other) semi-classical 
states. Usually this means that the semi-classical states are the states that we measure. 



the order of several times the Schwarzschild radius of the black hole. 
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4 Summary 



In this final section we give an executive summary of our main results. Consider the total 
quantum state of the black hole (with internal states | j ) ) after emission of quanta of late 
radiation (described by states | n)) 

l*> = -^E^U>l">l^>- (72) 
Here is the state of the early radiation. Here we assume that the late radiation started 



coming out at some time t = 0, and (72) is the state at some later time t = r. Since the 
early radiation is assumed to be dynamically decoupled, we are allowed to treat the states 
K&j) as static. 

We would like to investigate with which part of the wave function the radiation state 
is entangled, the black hole state or the early radiation state. In pQ it is argued that, 
since for an old black hole the Hilbert space of the early radiation is so much larger than 
that of the black hole itself, the radiation state must be maximally entangled with the 
environment state and has no entanglement with the black hole state |j )J^] This 
appears to exclude the existence of a smooth horizon. The question that we will first try 
to answer is: 

Under which conditions can an observer who has only access to the internal black hole 
states | j ) accurately determine in which state | the late radiation has been emitted? 

To focus the discussion, let us assume that an outside observer has detected the radiation 
in a particular state n). Thus instead of the full quantum state | ^ ) we now look at 

IW = 7rS^I J >l n >l $i >' (73) 

where we renormalized the state so that its norm is equal to one. Is there a measurement 
that acts on the states | j ) whose outcome is going to tell us about the quantum numbers 
n of the late radiation? 

The key point on which our paper is based is that the coefficients C^j are not just some 
arbitrary set of numbers: for given n they are transition amplitudes between black hole 
states, as a result of emitting the radiation in state | n) . Hence the C' ■ are determined by 
the microphysics, and represent the matrix elements of particular operators C n , that act on 



10 In section 2.4 we already established that this is not entirely true, even in the maximally entangled 
situation. The radiation state is always shares entanglement with the black hole. 
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the interior part of the black hole Hilbert space. We assume that the C^j are all statistically 
independent random matrices. Now pick some arbitrary state Iwi), with the same energy 
as | n) , and consider the operator U m = : C m l co d e : with matrix elements 

(U m ) jk = £ C ifU - trace (74) 

i 

Here w m = -|e~" m is the Boltzmann weight. The sum over % runs over some code subspace 
of the black hole Hilbert space of dimension N code < N. We will now show that this 
operator, which only acts on the black hole interior, is capable of recovering the information 
about the radiation state nV 



Applying this operator U m to ( 73 ) gives 



H m ||n» = lj2 C L ( U rn) jk \k)\n) | (75) 

Since the C^/s are ergodic matrices, whose coefficients all arbitrary complex phases, we 
can evaluate this sum with the help of standard statistical reasoning. We need to look for 



the resonant contributions. Inserting (74) into (75), we obtain an expression with a sum 
over two internal indices: j and % . The index j labels black hole states of mass M — E n and 
thus runs over N n indices and i labels the code subspace and thus runs over N code indices. 



We evaluate both sums by keeping only resonant terms, by using eqns like (38) and (65). 



We call the sum over j the 'good contraction' and the sum over i the 'bad contraction'. 



The second term in (74) is chosen such that it cancels the 'bad contraction'. The good 
contraction gives back the state |nj> when m = n, but with a projection on to the code 
subspace states 

n m ||n)) = <5 m JI code ||n), (76) 

where n co de = 5^„n n . So it looks like we have been able to establish that the black hole 
state is maximally entangled with the radiation mode \n}, provided that the black hole 
part of the initial state lies in the code subspace. So why can't we just take the code space 
equal to the complete interior Hilbert space and set iV co de = iV? 



There's the catch. Eqn (76) is true as a coarse grained equation. In computing the 
sums in ( 75 ) we ignored many off-diagonal terms, which individually are all very small. We 
should worry, however, about their collective effect. One place where this enters is in the 
verification that the operator Ii m really acts like a projection operator with Ii 2 m = Ii m . By 
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Figure 3: The QECC protects the coherence of a subspace % co de within the black hole Hilbert 
space %bh- It safeguards the smoothness of the horizon for every state in H co de- As the density 
matrix p spreads out, the code subspace needs to grow along with it. This degrades the fidelity 
of the code and spoils the semi-classical correspondence. The state decoheres into a sum of semi- 
classical components, each of which fits inside a smaller code space and has a smooth horizon. 



the same calculation method, we find 

n m = n m + E m , E m = codc i M -E m , (77) 

where lM-E m is the unit operator on the black hole Hilbert space with mass M — E m . The 
extra term E m the error of the recovery operation. So to have an accurate recovery without 
appreciable errors, we need that the code subspace is small. 

We conclude that the black hole is still maximally entangled with the radiation mode 
when the state of the black hole fits inside of some suitably chosen code subspace with 

n cod e|n)) = I n) and N code < N n . (78) 

In section 3.2 we show that in this case one can reconstruct the interior QFT observables in 
terms of the recovery operators. This is compatible with information theoretic constraints: 
the reconstruction works only as long as the entanglement entropy with the early radiation 
is sufficiently smaller than the entropy log N n of the space of the states | j ). 

However, we have also found that when the von Neumann entropy of the internal black 
hole state becomes too large, then the quantum recovery operation receives an error con- 



tribution of the form (77). The low energy effective field theory breaks down. This is 



consistent with the firewall argument. But we have learned more: with our methods we 
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have been able to qualify and quantify how the effective field theory description begins to 
break down as the entanglement with the external environment increases. 

What we have just described is a standard issue in quantum measurement theory. Sup- 
pose we think of the black hole as a measuring device. We can then ask under which 
circumstance can it measure the state of the radiation. Measurement in general happens 
when the measured object decoheres by entangling itself with the internal quantum states 
of the measuring device. The firewall argument invokes the situation where the complete 
state of a big measuring apparatus A (black hole) itself has been measured in the past 
by an even bigger system Z (early radiation). Monogamy of entanglement then seems to 
imply that the smaller apparatus A can no longer do any measurements on an even smaller 
system B (late radiation). This is not true, in general. 

It is useful to compare two situations: (a) the bigger system Z remains in sufficient 
contact with the smaller apparatus A, so that it keeps actively measuring the complete 
state of A, or (b) A was once maximally entangled with Z in the far past, but A and Z 
are dynamically decoupled at present. In case (a), A can indeed not do any measurements 
itself. This follows from the dynamics of dephasing: in other to measure some other small 
system, B, it needs to build phase coherence with B. But while A is being measured, 
all components of its state are randomized in ways that prevent its ability to build phase 
coherence with anything else. So in this case, entanglement is indeed monogamous. In case 
(b), however, the measuring apparatus A can in general still do measurements. Once it is 
dynamically decoupled from Z, it can start building new phase coherence with some other 
system B. So even if A started out in a maximally mixed state, over time it can build 
new entanglement via interactions with its environmen, and measure B. The black hole 
situation is more like case (b). 
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Appendix A: Some calculations 

In this Appendix we outline some of the calculations, the were left out of the main text in 



Chapter 3. We first do the calculation of the expectation value (54) of the product AB of 



an internal operator A given in (49) and an external operator B made out of b- modes. We 



first do an intermediate calculation 



6 (0| rfABU\0) b = b (0\ a (0\R 1 AR\0) a BU\o) t 
= b (0\ a (0\U^R^ABRU\0) a \0) b 



[0u\AB\0u) 1 



code- 



(79) 
(80) 
(81) 



Here we first use the definition (49 ), then move the B operators inwards (since they commute 



with the recovery operators), and finally use the result (42) for the action of RU on states 



of the form |z)|0) |0),. Here l co d e denotes the unit operator on "H co d e - Now we take the 



trace with the density matrix p(0) given in (47). Since p(0) is assumed to act within "H co de, 
we immediately find 



tv(p(r)AB) = J2Pik(i\(0b\rtABU\0 b )\k) 



i ,k 



= J2pik(t\(0u\AB0 u )\k) = (0 u \AB0 u ). 

i ,k 

This reproduces the expectation value in the local Minkowski vacuum. 



(82) 



(83) 



Next we will verify eqn (53). We will work in the limit that N CO( i e N. Using the 



definition (39), we have 



V w n w m _ 



34) 



We now compute 



RlR s R\R m = - 1 J2c n \i){j\CtC t \k){k\C, 



'm - - 
i ,fc 



>st 



(85) 
(86) 



3^ 



where we used eqn (38 



). In the same way, we derive (55) 



tr A (pRlR m ) = ^ Pli -^(i\C^C M \s)(s\ctC n \k) (87) 



\fw n w m 5 n t 5 ms , 



where we use (38) twice. 
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